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ABSTRACT 

We reanalyze the formation and evolution of galaxies in the hierarchical clustering scenario. Using a semi- 
analytic model (SAM) of galaxy formation described in this paper, which we hereafter call the Mitaka model, 
we extensively investigate the observed scaling relations of galaxies among photometric, kinematic, structural 
and chemical characteristics. In such a scenario, spheroidal galaxies are assumed to be formed by major 
merger and subsequent starburst, in contrast to the traditional scenario of monolithic cloud collapse. As a new 
ingredient of SAMs, we introduce the effects of dynamical response to supernova-induced gas removal on 
size and velocity dispersion, which play an important role on dwarf galaxy formation. In previous theoretical 
studies of dwarf galaxies based on the monolithic cloud collapse given by Yoshii & Arimoto and Dekel & 
Silk, the dynamical response was treated in the extremes of a purely baryonic cloud and a baryonic cloud fully 
supported by surrounding dark matter. To improve this simple treatment, in our previous paper, we formulated 
the dynamical response in more realistic, intermediate situations between the above extremes. While the effects 
of dynamical response depend on the mass fraction of removed gas from a galaxy, how much amount of the 
gas remains just after major merger depends on the star formation history. A variety of star formation histories 
are generated through the Monte Carlo realization of merging histories of dark halos, and it is found that our 
SAM naturally makes a wide variety of dwarf galaxies and their dispersed characteristics as observed. It is also 
found that our result strongly depends on the adopted redshift dependence of star formation timescale, because 
it determines the gas fraction in high-redshift galaxies for which major mergers frequently occur. We test four 
star formation models. The first model has a constant timescale of star formation independent of redshift. 
The last model has a timescale proportional to the dynamical timescale of the galactic disk. The other models 
have timescales intermediates of these two. The last model fails to reproduce observations, because it predicts 
only little amount of the leftover gas at major mergers, therefore giving too weak dynamical response on size 
and velocity dispersion of dwarf spheroidals. The models, having a constant timescale of star formation or a 
timescale very weakly dependent on redshift, associated with our SAM, succeed to reproduce most observations 
from giant to dwarf galaxies, except that the adopted strong supernova feedback in this paper does not fully 
explain the color-magnitude relation under the cluster environment and the Tully-Fisher relation. A direction 
of overcoming this remaining problem is also discussed. 

Subject headings: cosmology: theory - galaxies: dwarf - galaxies: evolution - galaxies: formation - large- 
scale structure of universe 



1. INTRODUCTION 

Dwarf galaxies give us many useful insights on the galaxy 
formation. They occupy a dominant fraction among galax- 
ies in number and have a wide variety of structures and for- 
mation histories. Furthermore it would be easier to under- 
stand how physical processes such as heating by supernova 
explosions affect their evolution, because of shallowness of 
their gravitational potential well, compared to massive galax- 
ies that have more complex formation histories and physical 
processes. Equally important is that massive galaxies have 
evolved via continuous mergers and accretion of less massive 
galaxies, according to the recent standard scenario of large- 
scale structure formation in the universe with the cold dark 
matter (CDM). This indicates that it is essential to understand 
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the formation of dwarf galaxies even in understanding the for- 
mation of massive galaxies, because they would be formed by 
mergers of their building blocks like dwarf galaxies. 

Before the CDM model gained its popularity, tra- 
ditional models of monolithic cloud collapse such 
as the galactic wind model for elliptical galaxies 
(Eggen, Lvnden-bell & Sandage 1962; Larson 1969; Ikeuchi 
1 19771 iSaitol 119791 I Yoshii & Said 119791; lArimoto & Yoshiil 
[19861 Il987t iKodama & Arimotdl 19971) and the infall m odel 
for spiral galaxies (Arimoto, Yoshii & Takahara 1991) had 
been widely used in analyses of galaxy evolution. These 
are simple but strong tools to explain various observations 
and then have contributed to construction of a basic picture 
of galaxy ev olution. Based on such traditional models, 
iDekel & Silk! (119861) systematically investigated many as- 
pects of dwarf galaxies. Especially they focused on the role 
of supernova (SN) feedback, that is, heating up and sweeping 
out of the galactic gas by multiple SN explosions, in the 
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processes of formation and evolution of dwarf galaxies. 
Since a large amount of heated gas is expelled by the galactic 
wind, the self-gravitating system exp ands as a result of 
dynamical respons e to the gas removal (Hills 1980; Mathieu 
19831 IVaderlfl986l) . They considered two limiting cases for 
the gas removal. One is a purely self-gravitating gas cloud, 
and another is a gas cloud embedded in a dominant dark 
halo. They showed that many observed scaling relations are 
well explained by their model. Combined with the evolu- 
tionary population synthesis code, lYoshii & Arimotol ( 119871) 
extended their analysis to estimating directly observable 
photometric properties, while only considering the purely 
self-gravitating gas cloud. 

Since the hierarchical clustering of dark halos predicted by 
the CDM model becomes a standard structure formation sce- 
nario, the galaxy formation scenario must be modified so as 
to be consistent with the hierarchical merging of dark halos. 
Explicitly taking into account the Monte Carlo realization 
of merging histories of dark halos based on the distribu- 
tion function of initial density fluctuations, the so-called 
semi-analytic models (S AMs) of galaxy formation have 
been developed ( e.g., iKa uffmann. White & Guiderdoni 
1991 lColeetal.1 11994 iBaueh Cole & FrenJd 119961 
NagaslrhTia"*Go"uda & SugiurdTT999 TlsT)mervi^ 
1999 1 : ISomerville. Primack & Faberl 120011 iNagashi ma et alJ 
20011 120021) . SAMs include several important physical 
processes such as star formation, supernova feedback, galaxy 
merger, population synthesis and so on. While SAMs well 
reproduce many observed properties of galaxies, most of 
analyses have been limited to those of massive galaxies. 
Thus we focus on the formation of dwarf galaxies in the 
framework of the SAMs in this paper. In the analyses of 
dwarf galaxies, the dynamical response to gas r emoval, which 
is known to play an important r ole as shown by D ekel & Silkl 
(1986) and Yoshii & Arimoto (1987), must be taken into 
account. While self-consistent equilibrium models taking 
int o account baryon an d da rk matter have been ana lyzed 
bv lYoshii & SaioT (fT987i) and ICiotti & Pellegrinil fl99l . we 
focus on the dynamical response of baryons within a dark 
matter halo. Since we have derived the mathematical formula 
of dynamical resp onse for the galaxies consi sting of baryon 
and dark matter in Nagashima & Yoshiil ( 120031) . it is possible 
to incorporate them into our SAM. 

The main purpo se of this paper is a reanalysis of 
iDekel & Silkl il986l) and lYoshii & Arimotol (119871) in the 
framework of SAMs. Thus we mainly focus on the forma- 
tion of elliptical galaxies, which are assumed to be formed 
by major merger and subsequent starburst. In order to do 
this, we construct the Mitaka model, which is a SAM in- 
cluding the effects of dynamical response. So far many 
scaling relations among photometric, structural, and kine- 
matical parameters of elliptical galaxies have been observed , 
such as the color-magnitude relation (e. g.. iBa urrJ 119591). 
the velocity dispersion-magnitude relation ( Faber & Jackson 



1976), and the surface brightness-size relation (Kormendv 
1977 1 iKodaira. Okamura & Watanabd U 9 8 3 ) . To understand 
how these relations are originated, especially on scales of 
dwarf galaxies, helps to clarify the processes of galaxy 
formation and evolution in the context of the cosmologi- 
cal structure formation scenario. Although some authors 
found by principal component analysis that elliptical galax- 
ies are distributed over the so-called fundamental plane in 
the three-dimensional space among photometric, structural 
and kinematical parameters ( Wata nabe. Kodaira & Okamural 



fT985llDresrier et alll987tlDiorgovski & Davisl 19871) . we fo- 
cus on the direct observables in this paper rather than their 
principal component projection. 

This paper is outlined as follows. In §2 we describe our 
SAM. In §3 we constrain model parameters in our SAM using 
local observations. In §4 we compare the theoretical predic- 
tions of SAM galaxies with various observations. In §5 we 
examine consistency check of our SAM with other observa- 
tions. In §6 we show the cosmic star formation history. In §7 
we discuss the nature of galaxies on the cooling diagram that 
has been traditionally used for understanding the formation of 
galaxies. In §8 we provide summary and conclusion. 

2. MODEL 

The galaxy formation scenario that we use is as follows. In 
the CDM universe, dark matter halos cluster gravitationally 
and merge in a manner that depends on the adopted power 
spectrum of the initial density fluctuations. In each of the 
merged dark halos, radiative gas cooling, star formation, and 
gas reheating by supernovae occur. The cooled dense gas and 
stars constitute galaxies. These galaxies sometimes merge to- 
gether in a common dark halo, and then more massive galax- 
ies form. Repeating these processes, galaxies form and evolve 
to the present epoch. 

Some ingredients of our SAM are revised. Modifications 
include the shape of mass function of dark halos, the star for- 
mation (SF) timescale, the merger timescale of galaxies tak- 
ing into account the tidal stripping of subhalos, and the dy- 
namical response to gas removal caused by starburst during 
major merger. The details are described below. 

2.1. Merging Histories of Dark Halos 

The merging histories of dark hal os are realized by 
a Mon te Carlo method proposed by Somerville & Kolatt 
( 1999) , based on the extended Press- S chechter (PS) formal - 
ism (lBondetan[T99lt iBoweil [T99ll lLacev & Colel fT993l) . 
This form alism is an extension of the Press-Schechter for- 
malism ( Press & Schecht erl 19741) . which gives the mass func- 
tion of dark halos, n(M), to estimate the mass function of 
progenitor halos with mass Mi at a redshift zo + Az of a 
single dark halo with mass Mo collapsing at a redshift zo, 
n(M\;zo + Az\Mo;zo)dMi. According to this mass function, 
a set of progenitors is realized. By repeating this, we obtain a 
merger tree. Realized trees are summed with a weight given 
by a mass function at output redshift. Dark halos with circular 
velocity V c i lc > Vi ow = 30km s -1 are regarded as isolated halos, 
otherwise as diffuse accreted matter. 

In our previous papers, we adopted the PS mass function to 
provide the weight for summing merger trees. Recent high- 
resolution Af-body simulations, however, sugge st that the PS 
mass function should be slightly cor rected (e.g.. | Jenkins et all 
120011) . According to the notation in lJenkins et alJ (120011) . the 
original PS mass function is written by 



2 5 C 

/(<7;PS)= \/--exp 
V 7r a 



2a 2 J 



(1) 



and the cumulative mass function n(M) is related with the 
above function by 

M dn(M) 



f(v,X) = 



po (finer 1 ' 



(2) 



where X specifies a model such as PS, po is the mean den- 
sity of the universe, a denotes the standard deviation of the 
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FIG. 1 . — Mass functions of dark halos. The s olid lin e represents the mass 
function given by Yahaai, Naaashima & Yoshii 1 2003, YNY), which is used 
in this paper. The dot-dashed line is the origin al PS mass functions, re spec- 
tively. The dotted line is the analytic result by Sheth & Tormen 1 1999, ST). 
The da shed line is the iV-body result over the resolved range by Jenk mset alj 
(2001, J). 



density fluctuation field, and S c is the critical density c ontrast 
for collapse assu ming spherically symmetric collapse ( Tomita 
U969tlGunn & Gotll972j). In this paper we use the fol l owing 
mass function given by Yahagi, Nagashima & Yoshii (2003, 
hereafter YNY) instead of the PS mass function, 



f(a; YNY) = A( 1 +x°)x D exp( 



-x 2 ), 



(3) 



where x = BS c /V2a, A = 2/[T(D/2) + T({C + D}/2)],B = 
0.893, C = 1.39 and D = 0.408. This is a fitting function that 
satisfies the normalization condition, that is, the integration 
over all rage of v is unity. These functions are plotted in Fig- 
ure[T]by the solid line (YNY) and the dot-das hed line (PS). We 
also sh ow other often used formula given by Sheth & Tormen 
( 1999), 
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and bv lJenkins 

/(cr;J) = 0.315exp(-| In cr" 1 +0.61 1 3 ' 8 ), 



(5) 



where A = 0.3222, a = 0.707 and p = 0.3 for the former (ST) 
and the latter (J) formula is valid over the range -1.2 < 
lncr-'<1.05. 

The number density of dark halos with large mass given 
by the above three functions (YNY, ST and J) and by the N- 
body simulation are similar to thos e given by analytic estima- 
tion (e.g., Yano, Nagashima & Gouda 1996; Monaco 1998; 
Nagashima 200 IJ) and are therefore trustworthy. With -body 
simulations, it is difficult to estimate the mass function at low 
mass because of the limited resolution and the uncertainty in 
the identification of dark halos. The very high resolution N- 
body simulatio ns, made using the adapt ive mesh refin ement 
method bv lYahagi & Yoshiil EoOll) and lYahagil (120021) . pre- 
dict a mass function similar to that of ST. This m ass function 
extend s to a mass 10 times smaller than that of J enkins et al.l 
d2001l) and has a slightly different slope. Thus we adopt the 
YNY mass function as above. A comparison between the 
YNY mass function and the jV-body simulation is dis cussed 
in a separate paper (lYahaei. N agashima & YoshIHl2003l) . In 



<j3] we show how the mass function of dark halos affects the 
luminosity function of galaxies. 

In short, we realize merger trees by using the extended PS 
formalism for dark halos whose mass function is given by the 
YNY at each output redshift. 

In this paper, we only consider a recent standard ACDM 
model, that is, fio = 0.3, fl\ = 0.7, h = 0.7 and a% = 0.9, where 
these parameters denote the mean density of the universe, the 
cosmological constant, the Hubble parameter and the nor- 
malization of the power spectrum of the initial density fluc- 
tuation fie l d. T he shape of the power spect rum given by 
Sueivama ( 1995), which is a modified one of Bardeen et alJ 
( 1986, hereafter BBKS) taking into account the baryonic ef- 
fects, is adopted. 

2.2. Tidal Stripping of Subhalos 

Most of recent high resolution A^-body simulations suggest 
that swallowed dark halos survive in their host halo as subha- 
los. Envelopes of those subhalos are stripped by tidal force 
from the host halo. We assume the radius of a tidally stripped 
subhalo r t by 



r t J" peri 



^"apo V^circ.h 



(6) 



where r pei i and r apo are the pericenter and apocenter for the 
orbit of the subhalo, respectively, and subscripts "h" and "s" 
indicate the host halo and subhalo, resp ectively. In this pa- 
per, a ratio of r nRr j/r ann = 0.2 i s assumed ( Ghi gna et al.ll 1 998t 
Okamoto & Habdl999l l2b00). Because a singular isothermal 
profile for subhalos is assumed, their mass decreases propor- 
tional to r t /r s . The mass of stripped subhalos is used in the 
estimation of dynamical response to gas removal on the size 
and ve locit y dispersion when satellite galaxies merge together 
(see§Q- 



2.3. Gas Cooling, Star Formation and Supernova Feedback 

The mean mass density in dark halos is assumed to be pro- 
portional to the cosmic mean density at the epoch of col- 
lapse using a spherically symmetric collapse model ( Tomita 
1969; Gunn & Gott 1972). Each collapsing dark halo con- 
tains baryonic matter with a mass fraction fib/fio, where 
fib is the baryon density parameter. We adopt a value of 
fl b h 2 = 0.02 that is recently suggested by the BOOMERANG 
Project mea suring the anisotropy o f the cosmic microwave 
background ( Netterfield et al. 2002). This value is an inter- 
mediate one between (0.64- 1.4) x 10~ 2 given by analysis of 
light element abundance produced by big bang nucleosynthe- 
sis ( Suzuki, Yoshii & Beersl l2000l) and (2.24 ±0.09) x 10" 2 
given by a nalysis of cosmic mi crowave background observed 
by WMAP ( Spergel et al. 2003). The baryonic matter consists 
of diffuse hot gas, dense cold gas, and stars. 

When a halo collapses, the hot gas is shock-heated to the 
virial temperature of the halo with an isothermal density pro- 
file. A part of the hot gas cools and accretes to the disk of 
a galaxy until subsequent collapse of dark halos containing 
this halo. The amount of the cold gas involved is calculated 
by using metallicity-depend ent cooling functions provided by 
Sutherland & Domta ( 1993). The difference of cooling rates 
between the primordial and metal-polluted gases is prominent 
at T ~ 10 6 K due to line-cooling of metals. Chemical enrich- 
ment in hot gas is consistently solved with star formation and 
SN feedback. The cooling is, however, very efficient in dark 
halos with a virial temperature of T ~ 10 6 K even in the case 
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of the primordial gas, so the metallicity dependence of cool- 
ing rate only slightly affects our results. In order to avoid the 
formation of unphysically large galaxies, the cooling process 
is applied only to halos with V c i rc < V cut =250 km s" 1 . This 
manipulation would be needed, because the simple isothermal 
distribution forms so-called "monster galaxi es" due to too ef- 
ficient cooling at the center of halos. While Cole et alJ 
adopted another isothermal distribution with ce ntral core in- 
stead of such a simple cutoff of the cooling and Benson et al.l 
(2003b) considered some additional mechanisms such as the 
heating of hot gas by SNe and by heat conduction from out- 
side as well as its removal by superwinds from halos, we take 
the above simple approach. The value of V cut is rather small 
compared with our previous paper and other SAMs. This is 
caused by our assumption that invisible stars have negligible 
fra ction , which is introduced to darken luminosity of galaxies 
(§ 32. 6> . This smaller value of V cut makes the color of large 
galaxies less red, which shows up on a bright p ortio n of the 
color-magnitude relation of elliptical galaxies (§ 35. 5> . 

Stars in disks are formed from the cold gas. The SF rate 
(SFR) M* is given by the cold gas mass M co id and a SF 
timescale r* as M* = M co id/r*. Now we consider two SF 
models. One is a constant star formation (CSF), in which t* 
is constant against redshift. Another is a dynamical star for- 
mation (DSF), in which t* is proportional to the dynamical 
timescale of the halo, which allows for the possibility that the 
SF efficiency is variable with redshift. We then express these 
SF timescales as 



r°[l + /3(y circ )] 
T*°[l+/3(V circ )] 



Tdynfe) 
Tdyn(0) 



(CSF), 
(l+zr(DSF), 



(7) 



where t° and a are free parameters, and indicates the ratio 
of the SF timescale to the reheating timescale by the SN feed- 
back defined by equation (|9) (see below). Pure DSF occurs 
when (7 = 0. Because the timescale of cold gas consumption 
is equal to r*/(l +/3— R), where R is the returned mass fraction 
from evolved stars {R = 0.25 in this paper), the mass fraction 
of cold gas in galaxies that is nearly constant against their 
magnitude is automatically adjusted by multiplyi ng (1 + 0). 
Hence the parameter a* originally introduced by ICole et all 
( 1994) is eliminated by introducing the factor (1 + 0). The 
parameter r° is so chosen as to match the mass fraction of 
cold gas with the observed fraction (see §3). T hereby the 
SF-rel ated parameters are constrained according to lCole et al.l 
(2000). Since not all of cold gas might be observed, the ob- 
served data give a lower limit to mass fraction of cold gas. 

Figure|2]shows the redshift dependence of SF timescale for 
the four SF models of CSF (dot-dashed line), DSF0 (a = 0; 
solid line), DSF1 (a = 0.5; dotted line) and DSF2 (a = 1; 
dashed line). Objects, which collapse at higher redshift, have 
higher density and therefore shorter dynamical timescale. Ev- 
idently, the DSF with smaller a gives more rapid conversion 
of cold gas into stars, compared with the CSF. The four SF 
models predict different mass fraction of cold gas at high red- 
shift, leading to quite different characteristics of dwarf galax- 
ies. This indicates, as will be clarified later, that the redshift 
dependence of SF timescale can be constrained particularly 
from observed structures and photometric properties of dwarf 
galaxies. 

Massive stars explode as Type II SNe and heat up the sur- 
rounding cold gas. This SN feedback reheats the cold gas at a 
rate of M re h ea t = A^coid / T re heat, where the timescale of reheating 
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FIG. 2. — Star formation timescales r* as a function of redshift z in a 
ACDM model. The timescale is normalized to unity at the present epoch. 
The dot-dashed line represent the constant SF model (CSF). The solid, dotted 
and dashed lines represent the variants of dynamical SF model with a — 
(DSF0), 0.5 (DSF1) and 1 (DSF2), respectively, according to the prescrip- 
tions in equationlTl 



is given by 



where 



^reheat — 



/3(Vcirc) 



Vhot 



(8) 



(9) 



The free parameters Vh ot and cthot are determined by matching 
the local luminosity function of galaxies with observations. 

With the above equations and parameters, we obtain the 
masses of hot gas, cold gas, and disk stars as a function of 
time or redshift. Chemical enrichment is also taken into ac- 
count, adopting the heavy-element yield of y = 2Zq, assuming 
the instantaneous recycling approximation with the returned 
mass fraction from evolved stars R = 0.25. All of newly pro- 
duced metals are released into cold gas, then by SN feedback, 
a part of them is expelled into hot gas. Some metals in hot gas 
are brought back to cold gas by subsequent cooling, and are 
accumulated in stars by their formation. 

2.4. Mergers of dark halos and galaxies 

When two or more progenitor halos have merged, the newly 
formed larger halo should contain at least two or more galax- 
ies which had originally resided in the individual progenitor 
halos. By definition, we identify the central galaxy in the 
new common halo with the central galaxy contained in the 
most massive one of the progenitor halos. Other galaxies are 
regarded as satellite galaxies. These satellites merge by ei- 
ther dynamical friction or random collision. The timescale 
of merging by dynamical friction is given by T m „ = /mrgTfnc, 
where Tf r ; c is given by Binnev & Tremaine ( 1987), which is 
estimated from masses of the new common halo and the 
tidally truncated subhalo. The parameter is set to 0.7 
in this paper. When the time elapsed after merging of a 
progenitor halo exceeds r lmg , the satellite galaxy is accreted 
to the central galaxy. On the other hand, the mean free 
ti mescale of random colli sion of satellite galaxies r co n is given 
bv lMakino & Hutl i ll 9971) . With a probability Af/r co u, where 
At is the time step corresponding to the redshift interval Az 
of merger tree of dark halos, a satell ite galaxy merges with 
anoth er satellite picked out randomly ( Somerville & Primack 
1999). 
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Consider the case when two galaxies of masses m\ and 
ma(> m\) merge together. If the mass ratio / = m\/m.2 is 
larger than a certain critical value of /bulge; we assume that 
a starburst occurs and that all of the cold gas turns into stars 
and hot gas, which fills the resulting halo, and all of the stars 
populate the bulge of a new galaxy. On the other hand, if 
/ < /bulge, no starburst occurs, and a smaller galaxy is sim- 
ply absorbed into the disk of a larger galaxy. Throughout this 
paper we use /b u i ge = 0.5, which gives a consistent morpho- 
logical fraction in galaxy number counts. 

2.5. Size of Galaxies and Dynamical Response to 
Starburst-induced Gas Removal 

We assume that the size of spiral galaxies is determined 
by a radius at which the gas is supported by rotation, un- 
der the conservation of specific angular momentum of hot 
gas that cools and contracts. We also assume that the ini- 
tial specific angular momentum of the gas is the same as that 
of the host dark halo. Acquisition of the angular momen- 
tum of dark halos is determi ned by tidal torques in the ini- 
tial density fluctuation field jWhitel|1984tlCatelan & Theunsl 
U996albHNagashima & Goudall998l) . The distribution of the 
dimensionless spin parameter A#, which is defined by A# = 
L\E\ X I 2 /GM 5 I 2 where L is the angular momentum and E is 
the binding energy, is well approxi mated by a log-normal dis- 
tribution <Mo. Mao & White! 19 98). 



p(\ H )d\ H = 



exp 



(In A H -In A) 



21 



2a\ 



din A w , (10) 



where A is the mean value of spin parameter and a\ is its log- 
arithmic variance. We adopt A = 0.03 and a\ = 0.5. When the 
specific angular momentum is conserved, the effective radius 
r e of a presently observed galaxy at z = is related to the initial 
radius Rj of the progenitor gas sphere via r g = ( 1 .68 /y/2)XnRi 
jFalil979tlFail & Efstathiou 1980; Fall 1983). The initial ra- 
dius Ri is set to be the smaller one between the virial radius of 
the host halo and the cooling radius. A disk of a galaxy grows 
due to cooling and accretion of hot gas from more distant en- 
velope of its host halo. In our model, when the estimated 
radius by the above equation becomes larger than that in the 
previous time-step, the radius grows to the new larger value 
in the next step. At that time, the disk rotation velocity Vd is 
set to be the circular velocity of its host dark halo. 

Size estimation of high-redshift spiral galaxies, however, 
carries uncertainties because of the larg e dispersion in their 
observed size distribution. For example, Simard et al. ( 1999) 
suggests only mild evolution of disk size against redshift, tak- 
ing into account the selection effects arising from the detec- 
tion threshold of surface brightness, although the above sim- 
ple model predicts disk size proportional to virial radius R vu - 
of host dark halos evolving as R vn oc 1/(1 +z) for fixed mass. 
Allowing for the possibility that the conservation of angular 
momentum is not complete, we generalize this size estimation 
by introducing a simple redshift dependence, 



1.68 

V2 



x H Ri(i+zy 



(id 



where p is a free parameter. We simply use p = 1 as a reference 
value in this paper. The effect of changing p emerges in the 
selection effects due to the cosmological dimming of surface 
brightness and in the dust extinction, because the dust column 
density also changes with galaxy size. This is discussed in |5] 



Size of early-type galaxies, which likely form from galaxy 
mergers, are primarily determined by the virial radius of the 
baryonic component. When a major merger of galaxies oc- 
curs, assuming the energy conservation, we estimate the ve- 
locity dispersion of the merged system. Now we assign the 
subscript to the merged galaxy, and subscripts 1 and 2 to 
the central and satellite galaxies, respectively, in the case of 
central-satellite merger, or to larger and smaller galaxies, re- 
spectively, in the case of satellite-satellite merger. Using the 
virial theorem, the total energy for each galaxy is 

Ei = - l -[M h V 2 + {M d +M coXi )V 2 l (12) 

where Mi, and Md are the masses of bulge and disk, respec- 
tively, and Vt and Vd are the velocity dispersion of bulge and 
the rotation velocity of disk, respectively. Assuming the virial 
equilibrium, the binding energy Eh between the progenitors 
just before the merger is given by 



E h : 



Then we obtain 



E t E 



1^2 



(M 2 /M l )E 1 + (M l /M 2 )E 2 ' 



E\ +E 2 + Eb = Eq. 



(13) 
(14) 



Just after the merger there is only the bulge component con- 
sisting of cold gas and stars in the merger remnant, whose 
velocity dispersion is directly estim ated from the a bove equa- 
tion. This procedure is similar to lCole et alJ (12000). although 
they argued it in terms of size estimation. Then the size of the 
system just after the merger is defined by 



GMi 



(15) 



where M,- = +M co id is the total baryonic mass of the merged 
system. 

Next, the cold gas turns into stars and hot gas. Newly 
formed stellar mass is nearly equal to M co id/(l + (3) and the 
rest of the cold gas is expelled from the merged system to the 
halo by SN feedback. The final mass after the mass loss, Mf, 
can be estimated from the known (3. Assuming the density 
distributions of baryonic and dark matters, the dynamical re- 
sponse on the structural parameters to the mass loss can be 
estimated. In this paper we adopt the Jaffe model ( Jaff ell983l) 
for baryonic matter and the singular isothermal sphere for 
dark matter, and assume slow (adiabatic) gas removal com- 
pared with dynamical timescale of the system. Defining the 
ratios of mass, size, density and velocity dispersion at final 
state relative to those at initial state by A4,R,Y and U, the 
response under the above assumption is approximately given 
by 



R 



l+D/2 
M+D/2 1 



U - 



l YRi + Df(z f )/2 
l+D/fe)/2 ' 



(16) 



(17) 



where M. = YR*, D = l/ytzf, y and z are the ratios of den- 
sity and size of baryonic matter to those of dark matter and, 
f(z) is a function defined in Appendix. The subscripts i and 
/ stand for the initial and final states in t he mass loss process. 
The de tails are shown in Appendix and Nagashi ma & Yoshul 
(2003). The parameter D indicates the contribution of dark 
matter to the gravitational potential felt by baryonic matter in 
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the central region of a halo. In actual calculation, we use a 
circular velocity at the center of dark halos, V ce nt defined be- 
low, to estimate D as 2V c 2 ent /V l 2 ; , instead of the ratios of size y,- 
and density z, which characterize the global property of dark 
halos. If there is negligible dark matter (D — > 0), the well- 
known result of adiabatic in variant (M* + M co \S)r emerges 
(e.g.. lYoshii & Arimotoll987h . In contrast, if there is negligi- 
ble baryonic matter (D — > oo), R and U become unity, that is, 
size and velocity dispersion do not change during mass loss. 
The effect of the dynamical response is the most prominent 
for dwarf galaxies of low circular velocity. 

Considering realistic situations, the baryonic matter often 
condenses in the central region and becomes denser than the 
average density in the halo. It is likely that the depth of 
gravitational potential well is changed when a part of bary- 
onic mass is removed due to SN feedback. To take into ac- 
count this process, we define a central circular velocity of 
dark halo Vcent- When a dark halo collapses without any pro- 
genitors, Vcent is set to Vcirc- After that, although the mass of 
the dark halo grows up by subsequent accretion and/or merg- 
ers, Kent remains constant or decreases by the dynamical re- 
sponse. When the mass is doubled, V ce nt is set to V C ir C at that 
time again. The dynamical response to mass loss from a cen- 
tral galaxy of a dark halo by SN feedback lowers V ce nt of the 
dark halo as follows: 

Vcent,/ _ M f /2 + M d (n/r d ) 

Vcent,/ ~ Mi/2+M d {ri/r d ) ' 

The change of V C ent in each time step is only a few per cent. 
Under these conditions the approximation of static gravita- 
tional potential of dark matter is valid during starburst. 

Once a dark halo falls into its host dark halo, it is treated as 
a subhalo. Because we assume that subhalos do not grow up 
in mass, the central circular velocity of the subhalos mono- 
tonically decreases. Thus this affects the dynamical response 
later when mergers between satellite galaxies occur. We ap- 
proximate that the resultant density distribution remains to be 
isothermal with Vcent at least within the galaxy size. 

2.6. Photometric Properties and Morphological 
Identification 

The above processes are repeated until the output redshift 
and then the SF history of each galaxy is obtained. For the 
purpose of comparison with observations, we use a stellar 
population synthesis approach, from which the luminosities 
and colors of model galaxies are calculated. Given the SFR 
as a function of time or redshift, the absolute luminosity and 
colors of individual g alaxies are calculated using a popula- 
tion synthesis code by Kodama & Arimoto ( 1997). The stel- 
lar metallicity grids in the code cover a range from Z* =0.0001 
to 0.05. Note that we now define the metallicity as the mass 
fraction of metals. The initial stellar mass function (IMF) that 
we adopt is the power-law IMF of Salpeter form, with lower 
and upper mass limits of 0.1M Q and 60M Q , respectively. 

In most of SAM analyses, it has been assumed that there 
is a substantial fr action of invisible stars such as brown 
dwarfs. iCole et alJ Ji99l introduced a parameter defined as 
T = (M] um + Mbd) / 'M\ um , where M\ um is the total mass of lumi- 
nous stars with mass larger than 0. 1M Q and Mbd is that of in- 
visible brown dwarfs. A range of T ~ 1-3 has be en assumed 
depend ing on ingredients of SAMs. For example. ICole et all 
(2000) assumed T = 3.07 in the case of J7b = 0.04. In this pa- 
per, however, we do not assume the existence of substantial 
fraction of invisible stars. If a large value of T is adopted, the 



mass-to-light ratio of galaxies is too hig h to agree with obser- 
vations. Thus we fix T = 1. iSomerville & Primackl (fl999h 
also adopted a small value of T = 1.25 in their ACDM.3 
model for fl h = 0.037 (in their notation / lum = 1 /T = 0.8). 

The optical depth of internal dust is consistently estimated 
by our SAM. We take the usual assumption that the abundance 
of dust is proportional to the metallicity of cold gas, and then 
the optical depth is proportional to the column density of met- 
als. Then the optical depth r is given by 

A^coldZcold ,-,„. 

rex , (19) 

n 

where r e is the effective radius of the galactic disk. There 
are large uncertainties in estimating the proportionality con- 
stant, but we ado pt about a factor o f two smaller value com- 
pared with that in C ole et all (|2000), otherwise it predicts too 
strong extinction to reproduce galaxy number counts, pre- 
sumably because our chemical yield is higher than theirs. 
Wavelength dependence of optical depth is assume d to be 
the sa me as the Galactic extinction curve given by ISeatonl 
( 1979) . Dust di stribution is simply assumed to be the slab 
dust llDisnev. Davies & Phillippsll9 89). according to our pre- 
vious papers. We found that the resu ltant mean extinct ion for 
spiral galaxies is close to a model bv ICole et alJ {2000). 

We classify galaxies into different morphological types ac- 
cording to the B-band bulge-to-disk luminosity ratio B/D. In 
this paper, following Sim ien & de Vaucouleursl lf l9 86). galax- 
ies with B/D > 1.52, 0.68 < B/D < 1.52, and B/D < 0.68 
are classified as elliptical, lenticular, and spiral ga laxies, 
respectively. [K auffmann. White & Guiderdoni ( 19921) and 
Baug hTCole & Frenkl (119 96) showed that this method of type 
classification well reproduces the observed type mix. 

3. PARAMETER SETTINGS 

As already mentioned, we adopt a standard ACDM model. 
The cosmological parameters are 51o = 0.3, J7a = 0.7, h = 0.7 
and erg = 0.9. The baryon density parameter £lb = 0.02/z -2 is 
used. 

The astrophysical parameters are constrained from lo- 
cal observations , according to the procedure discussed in 
Nagashim a'et alJ (1200 ll 120021) . The adopted values are 
slightly different from those in our previous papers. This is 
mainly c ause d by adopting the differ ent m ass function of dark 
halos (§ 32. 1> and by fixing T = 1 (§ 32. 6> . Values of these pa- 
rameters are tabulated in Tables^andlJ] 

First, the SN feedback-related parameters (Vh t, ckta) and 
merger-related parameter (/nirg) are almost uniquely deter- 
mined if their values are so chosen as to reproduce the local 
luminosity function. Figure [3] shows theoretical results for 
CSF (solid line), DSF2 (dashed line), DSF1 (dot-dashed line) 
and DSF0 (dotted line). As in our previous papers, the SF 
timescale affects the local luminosity function only slightly. 
Symbols with errorbars represent observational results from 
the fi-band r edshift surveys, such as Automatic Plate Ma- 
chine (APM: lLovedav et al.l ll992). ESO Slice Project (ESP; 
IZucca et alJll997l). Durham/United Kingdom Schmidt Tele- 
scope (UKSTJRatcliffe etal] GUI) and Two-Degree Field 
(2dF; iFolkes et al.||1999 |). and from the /T-band redshift sur- 
veys given by Szokolv et al. ( 1998), Two Micron All Sky Sur- 
vey (2MASS; Koch aneketalJl2001l) and 2dF combined with 
2MASS (Cole etalJ2001l) . 

In Figure|3]we see how the mass function of dark halos af- 
fects the luminosity function. In the same figure we show 
the same model as CSF but for the PS mass function rep- 
resented by the thin solid lines. As the PS mass function 
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TABLE 1 

Standard Settings of Astrophysical Parameters in the Mitaka Model 



Parameter Value Annotation Observation 



^ lo1 150km s \ SU p ernova feedback-related (§§2.3) luminosity functions (Figure 3) 

c*hot 4 J 

Vnu 250 km s cooling cut-off (§§2.3) luminosity functions (Figure 3) 

V\ ow 30 km s minimum circular velocity of dark halos (§§2.1) — 

y 2Zq heavy-element yield (§§2.3) metallicity distribution (Figure 18) 

/bulge 0.5 major/minor merger criterion (§§2.4) morphological counts 

/mrg 0.7 coefficient of dynamical friction timescale (§§2.4) luminosity functions (Figure 3) 

<r 5^ spin parameter distribution (§§2.5) disk size (Figure 5) 

p 1 redshift dependence of disk size (§§2.5) faint galaxy number counts (Figures 20 and 21) 

T 1 fraction of invisible stellar mass (§§2.6) mass-to-light ratio (Figure 17) 



NOTE. — Cosmological parameters are: (Q , S^A . h, ct 8 , O b ) = (0.3, 0.7, 0.7,0.9, 0.02ft 2 ). 



TABLE 2 
Star 

Formation-Model 
Dependent 
Parameters 



Model 


r° (Gyr) 


<7 


CSF 


1.3 




DSF0 


1.7 





DSF1 


1.3 


0.5 


DSF2 


1.5 


1 



NOTE. — is determined 
by matching the mass frac- 
tion of cold gas with the ob- 
served fraction in spiral galax- 
ies reQuationlTland Figure 4]. 



predicts more dark halos and hence more galaxies than the 
YNY mass function except for the largest mass scale. Thus, 
we need Vh ot = 150 km s , less than 280km s" 1 in the previ- 
ous model, to weaken SN feedback which mainly determines 
the scale of the exponential cu t-off of the luminosity function 
(Nagashim a et aljl2001l I2002T) . We also investigated the ef- 
fects of power spectrum of density fluctuations and confirmed 
that neither a erg = 1 model nor a model with BBKS power 
spectrum without the baryonic effects is significantly differ- 
ent from our reference model. In we discuss the slight 
effect that <j% has on high-redshift galaxies. 

The effects of photoionization on the luminos- 
ity function have been discussed (iChiba & Nathl 
19941 iNagashima. Gouda & Sugiural Il999t 

Bullock. Kravtsov & Weinberg 2000; Nagashima & Gouda 
200 It ISomervilld l2002t iTulivetall 120021 iBensonetalJ 



2002a b, 2003a). There are two ways to suppress the gas 
cooling by photoionization. One effect is that the ultraviolet 
(UV) background from quasars and/or young stars prevents 
hot gas at the outer envelope from cooling. Another is that 
the Jeans mass becomes larger after cosmic reionization, 
which means larger Vi ow . It has been found that both pro- 
cesses lower the faint-end of the luminosity function. For 
example, if we use a large value of Vi ow ~ 70 km s" 1 for our 
model, which corresponds to increasing the Jeans mass, the 
number of galaxies at the faint-end of the resultant luminosity 
function decreases by about a factor of two, less than that 
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FIG. 3. — Local luminosity functions in the (a) B band and (b) K band. 
The thick solid, dashed, dot-dashed, and dotted lines represent the models of 
CSF, DSF2, DSF1, and DSF0, respectively, based on the YNY mass function. 
Only for the purpose of comparison here, shown by the thin solid line is the 
CSF model based on the original PS mass function. Symbols with error bars 
in (a) indicate the observational data from APM (Loveday et al. 1992, filled 
squares), ESP (Zucca et al. 1997, open triangles), Durham/UKST (Ratcliffe 
et al. 1998, open diamonds), and 2dF (Folkes et al. 1999, open circles). 
Symbols in (b) indicate the data from Szokoly et al. (1998, open squares), 
2MASS (Kochanek et al. 2000, open stars), and 2dF combined with 2MASS 
(Cole et al. 2000, crosses). 



given by the APM survey at -20 < M B -5\og(h) < -17 
(see Figure [3}- In that case we need weaker SN feedback, 
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FIG. 4. — Cold gas mass relative to B-band luminosity of spiral galaxies. 
The solid line represents the CSF model. Other SF models, not shown, give 
almost the same results with only 0.1 dex difference. The open squares in- 
dicate the observational data for atomic neutral hydrogen taken from Hucht- 
meier & Richter (1988). Since the cold gas in the model consists of all species 
of elements, its mass is multiplied by 0.75, i.e., Mm = 0.75-M co id, which cor- 
responds to a total fraction of atomic and molecular hydrogen. Therefore, 
in comparison with the theoretical result, the observational data should be 
regarded as lower limits. 



that is, smaller VW and/or smaller oth ot . We found, how- 
ever, that such a large value of V] ow does not make dwarf 
spheroidals at M B -5\og(h) > — 12 for Vi ow = 70 km s" 1 and at 
M B -5\og(h) > -10 for Vi ow = 50 km s" 1 . In addition, giant 
galaxies obtain more gas that has not cooled in smaller halos 
Vcirc < Viow, causing the bright-end of the luminosity function 
to shift brighter. Note that this is not effective for Vi ow < 40 
km s" 1 . Keeping in mind that such effects might affect our 
analysis through the determination of SN feedback-related 
parameters, we use Vi ow = 30 km s" 1 in this paper. 

Next, the SFR-related parameter (r°) is determined by us- 
ing the mass fraction of cold gas in spiral galaxies. The gas 
fraction depends on both the SN feedback-related and SFR- 
related parameters. The former parameters determine the gas 
fraction expelled from galaxies and the latter the gas fraction 
that is converted into stars. Therefore, in advance of deter- 
mining the SFR-related parameters, the SN feedback-related 
parameters must be determined by matching the local lumi- 
nosity function. 

Figure|4]shows the ratio of cold gas mass relative to B-band 
luminosity of spiral galaxies as a function of their luminos- 
ity. Theoretical result is shown only for the CSF model by 
the solid line. Other SF models of DSF2, DSF1 and DSF0 
provide almost the same results, and their differences from 
the CSF are only about 0.1 dex. We here assume that 75% 
of the cold gas in the models is compri sed of hydrogen, i.e., 
Myr_= 0.75M co id. H I data, taken from Huchtmeier & Richter 
( 1988), are shown by open squares with errorbars. Since their 
data do not include the fraction of H2 molecules, they should 
be regarded as providing a lower limit to the mass fraction of 
cold gas. Adopted values of the parameters in the SF models 
are tabulated in Table [2] Slight difference in the values of t° 
stems from the different redshift-dependence of SF timescale. 

Figure |5] shows the effective disk radii of local spiral 
galaxies as a function of their luminosity only for the CSF 
model. Other SF models also provide almost the same re- 
sults as CSF and are not shown. Thus, with the use of 
(A,ita)=(0.03, 0.5), all the SF models well reproduce the 
observe d disk size-m agnitude relation (thin solid line) com- 
piled bv lTotani & Yosbti J2000l) based on the data taken from 
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FIG. 5 . — Disk size of spiral galaxies. The thick solid line connecting open 
circles shows the theoretical result for effective radii of spiral galaxies in the 
CSF model. Errorbars on this line is the la scatter in predicted sizes. Other 
SF models, not shown, give almost the same results with less than 0.1 dex 
differe nce. The thin line indicates the mean relation given by Totani & Yoshii 
1 2000 ) for spiral galaxies, based on the observational data (crosses) taken 
from Impey et al. (1996). 



Imp ev et alJ §996), while showing a slightly steeper slope 
than the observed one. 

4. RESULTS 

4.1. Relationships between Size and Magnitude 

First, we examine the size distribution of elliptical galax- 
ies, because it is directly affected by the dynamical response 
to gas removal. Figure [6] shows the contour distribution in 
the effective radius versus absolute B-magnitude diagram. 
The levels of contours in order from outside to inside are 
0.02, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8 and 0.9 times 
the largest number of galaxies in grids. This way of de- 
scription applies to all the contour distributions shown be- 
low. The four panels show the results for the models of 
CSF, DSF2, DSF1 and DSF0, respectively, as indicated. Ef- 
fective radius r e of model galaxies are defined by 0.744r£, 
where rj is the three-dimensional hal f-mass radius, assum - 
ing the de Vaucouleurs profile (Nagas hima & Yoshiil l2003). 
The fill ed trian gles and squares represe nt the data t aken f rom 
iBender. Burstein & Faberl (1 1 9921 fl993h and lMateoH 19981) . re- 
spectively. When some of dwarf spheroidals are listed in 
both their tab les, we use the data in Bender. Burst ein & Faberl 
i 19921 119931) . Observed elliptical galaxies have two se- 
quences of dispersed dwarf galaxies with low surface bright- 
ness and compact dwarf galaxies with high surface brightness. 

As shown in Figure |2] the CSF model has the longest SF 
timescale at high redshift among the four SF models. This 
means that when a major merger occurs in the CSF model, a 
significant fraction of gas is removed, so that the dynamical 
response of dwarf galaxies is expected to be the largest. Fig- 
ure|6]clearly shows that the CSF model (panel a) predict many 
dwarf galaxies with very large size of r e ~ lO-lOVkpcfor 
M B -5 log/t > -10. In the range of -10 > M B -5\ogh > -15, 
the peak of distribution is on the sequence of dispersed dwarf 
ellipticals. On the other hand, the DSF0 model (panel d) has 
the shortest SF timescale at high redshift, so that the effects 
of gas removal are not significant. Most of galaxies are dis- 
tributed along a single power-law sequence, corresponding to 
the compact dwarf galaxies. The slop e is determined by ahot 
as shown in Nagashima et al. (2001), in which they simply 
determined the size of elliptical galaxies by GMb/2V£. The 
behavior of SF timescale indicates that the models of DSF1 
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: (a) CSF 



:: (b) DSF2 
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FIG. 6. — Distribution of effective radius for elliptical galaxies in the mod- 
els of (a) CSF, (b) DSF2, (c) DSF1, and (d) DSFO, where the effects of dy- 
namical response to starburst-induced gas removal are taken into account. 
Based on the calculated number of such galaxies per magnitude per loga- 
rithm of kpc, the levels of contours shown in order from outer to inner are 
0.02, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, and 0.9 times the largest number of 
galaxies in grids, respectively Symb ols indicate th e observationa l data given 
by Bender, Burstein & Fabei 1 1992) (triangles) and Mateo 1 1998) (squares). 
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FIG. 7. — Distribution of effective radius for elliptical galaxies. Same as 
Figurelol but for the SF models without the effects of dynamical response to 
starburst-induced gas removal. 



and DSF2 are intermediate between CSF and DSFO. This 
suggests that the SF timescale should be constant indepen- 
dent of redshift or very weakly dependent on redshift, but not 
proportional to the dynamical timescale of galaxies. 

It should be noted that the gas removal is not a unique mech- 
anism to make large ellipticals. If a galaxy just after major 
merger keeps a large fraction of gas, its size becomes larger 
than that of gas-poor system with the same stellar mass be- 
cause the size is proportional to the total mass of stars and 
cold gas. This can be seen directly if we do not take into ac- 
count the effects of dynamical response to gas removal. Fig- 
ure shows the same result as Figure [6] except without the 
effects of dynamical response. The overall shapes are similar, 
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FIG. 8. — Distribution of B-band surface brightness within the effective 
radius for elliptical galaxies in the models of (a) CSF, (b) DSF2, (c) DSFl, 
and (d) DSFO, where the effects of dynamical response to starburst-induced 
gas removal are taken into account. Contours and symbols have the same 
meanings as in Figurelol 



independent of SF timescale, because the effects of dynamical 
response are not taken into account. However the peak loca- 
tions in Figure [7] are different from those in Figure [6] This 
is because the mass of cold gas just after major merger deter- 
mines its size as mentioned above and therefore the size de- 
pends on the SF timescale. In other words, the CSF model has 
a larger fraction of cold gas and hence a larger galactic mass at 
high redshift, compared with the DSF models. This situation 
without the effects of dynami cal response co rresponds to the 
case of dominant dark halo in Dekel & Silk ( 1986), in which 
the size and velocity dispersion do not change during the gas 
removal. We therefore conclude that in realistic situations the 
dark matter is not always dominant in the gravitational poten- 
tial. Note that while the r e -Ms relation is reproduced even 
without the effects of the dynamical response, as will shown 
in the next subsection, the Faber-Jackson relation is not repro- 
duced unless the effects are properly taken into account (see 
FigureslT21andrnV 

In addition to size and magnitude, surface brightness is 
also an important observable quantity. Two sequences of 
dispersed and compact dwarf galaxies are prominent in the 
surface brightness versus absolute fi-magnitude diagram, as 
shown in Figure [8] Surface brightness of model galaxies is 
defined as the average brightness in the area encircled by ef- 
fective radius. For giant ellipticals, surface brightness is pre- 
dicted to become brighter towards brighter magnitude with 
shallow slope, converging to ii e $ ~ 22, as observed. For 
Mg-5\ogh > -18, surface brightness widely spreads. This 
magnitude corresponds to Vh ot , that is, the magnitude at which 
j3 becomes larger than unity. The CSF and DSF2 models 
make the widely spread distribution which reproduces two se- 
quences simultaneously and thus are likely, while the DSFO 
makes one sequence of compact ellipticals only. 

In all the SF models there are many galaxies with very 
low surface brightness for Mb- 51og/i > — 15. Actually such 
galaxies cannot be detected unless the detection threshold of 
surface brightness is faint enough in galaxy survey observa- 
tions. With a cutoff of ii c ^b = 26.5 taken into account in both 
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FIG. 9. — Distribution of B-band surface brightness within the effective radius for elliptical galaxies. Same as Figure|S] but for only galaxies of high surface 
brightness with /i e g < 26.5 in the models as well as the data. 



the model and observed data, we show the surface brightness 
versus absolute B-magnitude diagram in Figure [9] and the ef- 
fective radius versus absolute B-magnitude diagram in Figure 
fTol For the models of CSF and DSF2 the resultant distribu- 
tions in both of these diagrams are consistent with the data. 

As another combination of observable quantities, the sur- 
face brightness vers us size rela tion is often referred to as the 
Kormendy relation ( Kormendvl l 1 9771) . Figure fTTI shows the 
theoretical distribution of elliptical galaxies with /vb < 26.5 
in the surface brightness versus effective radius diagram. Ex- 
tended distribution toward low surface brightness is clearly 
seen for the models of CSF and DSF2. 

4.2. Faber-Jackson Relation 

Velocity dispersion of elliptical galaxies is also an indepen- 
dent dynamical observable. If the dark matter dominates the 
baryonic matter, the velocity dispersion of galaxies reflects 
that of their host dark halos. On the other hand, in the case of 
negligible dark matter, the velocity dispersion is substantially 
affected by dynamical response to gas removal in proportion 



to the mass fraction of removed gas. Thus the velocity disper- 
sion versus absolute magnitude relation for elliptical galax- 
ies, often called Faber-Jackson (FJ) relatio n, provides an- 
other strong constraint on galaxy formation ( Fabe r & Jackson! 
1976). 

Figure [2] shows the velocity dispersion versus absolute 
B-magnitude diagram for the models of CSF, DSF2, DSF1 
and DSFO. Velocity dispersion of model galaxies is assumed 
to be isotropic and is converted to one-dimensional central 
dispersion by err/ID) = Vb/\/3 after increased to the central 
value by a factor of V2 according to the de Vaucouleurs-like 
profile. For the DSFO model, galaxies with velocity disper- 
sion less than 10km s" 1 is scarcely distributed. Note that 
the cutoff circular velocity Vi ow , above which dark halos are 
identified as isolated objects, is 30km s -1 nearly correspond- 
ing to the Jeans scale for collapse after cosmic reionization 
(Ostri ker& GnedinTT9 96). This makes the sequence that con- 
verges to CT (1D) = 30/^ ~ 17 km s" 1 . For the DSFO model, 
we see only such sequence because the smallest fraction of 
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FIG. 10. — Distribution of effective radius for elliptical galaxies. Same as Figure[S| but for only galaxies of high surface brightness with fi e g < 26.5 in the 
models as well as the data. 



gas at major merger in this model gives very weak effects of 
dynamical response to gas removal. It should be noted that 
this tendency hardly depends on Vi ow because the SN feed- 
back is very efficient for dwarf galaxies. In contrast, other SF 
models can reproduce the low velocity dispersion observed 
for local dwarf spheroidals for M B - 5 log h > - 1 5 . 

It is possible that dwarf galaxies with very low velocity dis- 
persion would remain undetected from actual observations, 
because such galaxies are expected to have very low sur- 
face brightness as a result of dynamical expansion by gas re- 
moval. Figure fT~3l therefore shows the FJ relation after ex- 
cluding galaxies of lower surface brightness with /j, e> B > 26.5. 
Most of dwarf galaxies of low circular velocity are excluded 
from the CSF model and their distribution is in apparent dis- 
agreement with the data. This situation is a little improved for 
the models of DSF1 and DSF2. Figure [T3l clearly indicates 
the difficulty of making dwarf galaxies of low circular veloc- 
ity and high surface brightness. In other words, the FJ rela- 
tion strongly constrains the SF timescale requiring the mild 
or negligible redshift dependence, when compared with the 



dynamical timescale. 

Here we examine the effects of the dynamical response on 
velocity dispersion. Figure PHI shows the velocity dispersion 
versus absolute B-magnitude diagram without the effects of 
dynamical response as in Figure Q Evidently, all four SF 
models cannot reproduce the dwarf galaxies of low circular 
velocity and give almost the same distribution on this dia- 
gram independent of SF timescale. This hilights the impor- 
tance of dynamical respons e for th e formation of dwarf galax- 
ies. Kauffmann & Chariot ( 1998) predicted the FJ relation in 
their SAM with no dynamical response taken into account and 
claimed to find good agreement with observations. Their re- 
sult is, however, limited only to Mb < -18, where such dy- 
namical response has no significant effect. Thus our SAM 
analysis is the first that has reproduced the observed velocity 
dispersion even for dwarf galaxies. 

4.3. Mass-to-Light Ratio 

The mass, which is used to estimate the mass-to-light ra- 
tio, is the dynamical mass of Md yn oc r^cr 2 /G. In our SAM, 
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FIG. 1 1. — Kormendy relation for elliptical galaxies with fi e s < 26.5 in the models of (a) CSF, (b) DSF2, (c) DSF1, and (d) DSFO, where the effects of 
dynamical response to starburst-induced gas removal are taken into account. Contours and symbols have the same meanings as in Figure l9l 



the size and the velocity dispersion a are determined by 
the amount of baryonic matter that escapes out of the system 
against the underlying dark matter potential. In this paper we 
define the dynamical mass by M dyn = 2r b Vg/G. Figure El 
shows the mass-to-light ratio of elliptical galaxies as a func- 
tion of absolute B-magnitude for the four SF models. 

Our SAM well reproduces the observed mass-to-light ratio 
for giant ellipticals, and at least qualitatively the mass-to-light 
ratio for dwarf spheroidals that increases towards faint mag- 
nitude. As for giant ellipticals, observed data still have a large 
scatter. This might be caused by an uncertainty in estimating 
dynamical mass owing to anisotropic kinematics and rotation. 

Such trend of the mass-to-light ratio for dwarf spheroidals 
is not caused by the dynamical response alone. Figure El 
shows the mass-to-light ratio with no dynamical response 
taken into account, which turns out to be very similar to or 
slightly higher than that with dynamical response in Figure 
IT31 Therefore, considerably high values up to M /Lb ~ 10 3 
are caused by the dominant dark halo in which the size and 
velocity dispersion are kept unchanged during the starburst. 



Since the luminosity Lb is proportional to the final mass My, it 
follows thatM/L B oc M,-/My oc (3 in the limit of j3 > 1, where 
j3 measures the strength of the SN feedback in equation 0. 
Using a relation M oc Vj„ from th e spherical collapse model 
( Tom itall969tlGunn & Gotil972l we obtain 

— oc (3 oc <j- Qta oc M - ahJ(3+ah ° l) . (20) 
Lb 

If we adopt a hot = 4, then M/L B oc M~ Ap oc 10 8M "/ 35 . The 
slope of this relation explains the result in Figure El There- 
fore, the small er mass fraction of d a rk matter lowers the mass - 
to-light ratio JDekel & SiUdll986t lYoshii & Arimotolll987l) . 
In apparent contrast to this, the DSFO model scarcely forms 
galaxies of such high mass-to-light ratio. This is because 
the size is determined by the baryonic mass, r e oc GMb/V^, 
and because the removed gas is negligible due to the short SF 
timescale at high redshift. 

By excluding galaxies of lower surface brightness with 
He,B > 26.5, we obtain gross agreement with the observed 
mass-to-light ratio in FigurefTTl 
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FIG. 12. — Faber-Jackson relation for elliptical galaxies in the models 
of (a) CSF, (b) DSF2, (c) DSF1, and (d) DSFO, where the effects of dy- 
namical response to starburst-induced gas removal are taken into account. 
One-dimensional central velocity dispersion cto(ID) for model galaxies is es- 
timated from Vt I v3 after increased to the central value by a factor of v2 
according to the de Vaucouleurs-like profile. Contours and symbols have the 
same meanings as in Figurelol 



4.4. Metallicity-Magnitude Relation 

Figure ^] shows the mean stellar metallicity as a function 
of absolute fi-magnitude for elliptical galaxies. The verti- 
cal axis represents the logarithmic iron abundance [Fe/H]. 
Shown by symbols are the data compiled by M ated Jl99Sh 
and those translated from Mg? index for those taken by 
iBender. Burstein & Faberl i ll 9931) . whereas by contours the 
Lg-weighted average of logarithmic metal abundances of stars 
(log(Z*/Z Q ))£ B for the theoretical prediction. Since the theo- 
retical results are almost independent of whether galaxies of 
low surface brightness are excluded, we only show the metal- 
licity distribution excluding such galaxies with fi e g > 26.5. 
We do not show the metallicity distribution for spirals, which 
is similar to that for ellipticals. 

All the SF models well reproduce the observed tight rela- 
tion between metallicity and absolute magnitude. This means 
that such relation is not affected by the redshift dependence 
of SF timescale, because we adopt the SN feedback in the 
same way irrespective of either continuous or burst forma- 
tion of stars in galaxies. With strong SN feedback with 
«hot — 4, our SAM succeeds to explain the observed low 
metallicity [Fe/H] ~ -2 at Mb - 5 log h ~ - 1 0, in spite of using 
a rather high value of metallicity yield y = 2Z Q . Note that if 
ahot = 2, we obtain about a factor of three higher metallicity 
at that magnitude . This is consistent with the conclusion by 
iCole et aIT ( l2000l) that with ahot = 2 the metallicities of their 
dwarf galaxies are systematically higher than the data, while 
such predicted values reside in a range of large observed scat- 
ter. 

5. CONSISTENCY WITH OTHER OBSERVATIONS 

In this section, we show the whole aspects of our SAM to 
check the consistency with other local and high-redshift ob- 
servations. This will clarify the limitations of present SAM 
analyses and physical processes for which further investiga- 
tion is required. In the first three subsections (§§5.1-5.3), 



high-redshift properties of model galaxies are examined by 
comp aring the observed ga laxies in the Hubble Deep Field 
(HDF.IWilliams et alll996l) and the Subaru Deep Field (SDF, 
Maih araetalJl2001l) . Then, the Tully-Fisher relation (TFR) 
for local spiral galaxies (§5.4) and the color-magnitude rela- 
tion (CMR) for cluster elliptical galaxies (§5.5) are discussed. 

5.1. Faint Galaxy Number Counts 

Counting galaxies as a function of apparent magnitude is 
one of the most important observable quantities for constrain- 
ing th e geometry of the univers e and the evolution of galaxies 
(e.g.. lYoshii & Takahardfl9 881 We already showed that our 
SAM can simultane ously reproduce the UV /optical galaxy 
counts in the HDF (Nagashima et al. 2001) and the near- 
infrared galaxy counts in the SDF llNag ashima et al J 120021) . 
These previous works demonstrated that inclusion of the se- 
lection effects arising from the cosmological dimming of sur- 
face b rightness of high-redshift galaxies derived by lYoshiil 
i ll 9931) is essential, because the number count of galaxies 
is obtained by summing up the product of luminosity func- 
tion and cosmological volume element to the accessible max- 
imal redshift above which galaxies have the surface bright- 
ness fainter than the threshold and thus are not detected. De- 
ta ils of the meth o d to estimate such selec ti on effects are found 
in lYoshiil J1993I) iTotani & Yoshiil J2000l) . lTotani et all J2001I) 
and lNagashima et alJ l l200lU2002l) . 

In this paper the YNY mass function of dark halos (§ 32. 1> 
is used instead of the PS, and the dynamical response to 
gas removal is a novel ingredient in the analysis. With 
these modifications we reexamine the galaxy counts again. 
Figure [191 shows the UV/optical galaxy counts in the HDF 
(£7300,^450,^606 and in the AB system). For reference, 
other ground-based observed counts are also plotted. The 
thick solid and dashed lines show the predictions of CSF 
and DSF2, re spectively, taking into account the absorption by 
internal dust (Totani & Yoshii 2000) and by intervening H I 
clouds ( Yoshii & Peterson 1994), and the selection effects due 
to cosmological dimming of surface brightness (Yoshii 1993). 
We adopt p = 1 in equation ill It . which determines the disk 
size at high redshift. The thin lines represent the same mod- 
els, except without the selection effects. We do not show the 
models of DSF1 and DSFO, because they are much the same 
as CSF and DSF2 with only small difference comparable to 
observational errors. These models only slightly underpredict 
th e tAoo count and o verpredict the 7si4 count, as pointed out 
byEagashima et alJ J200 lh . 

Figure|20|shows the near-infrared galaxy counts in the SDF 
(K 1 ). Like the HDF counts, the models of CSF and DSF2 well 
reproduce the SDF counts. For the purpose of comparison, 
the DSFO model is also shown. Because of the highest SFR at 
high redshift, the stellar mass in this model increases rapidly, 
making too many faint galaxies to reproduce the observation. 
Therefore, the models of CSF and DSF2 give the predictions 
that agree with the HDF and SDF counts simultaneously. 

5.2. Redshift Distribution 

Figure 1271 shows the redshift distributions of /^-selected 
galaxies in the HDF with 22 < hu < 24 (panel a), 24 < 
/gi4 < 26 (panel b) and 26 < 4i4 < 28 (panel c). The num- 
ber of model galaxies in each panel is calculated over the 
same celestial area as the HDF. The thick solid, dashed and 
dot-dashed lines show the results of CSF, CSF with erg = 1 
and DSFO, respectively. The histogram in e ach panel is the 
observed photometric redshift distribution bv lFurusawa et alJ 
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FIG. 13. — Faber-Jackson relation for elliptical galaxies. Same as Figure [T2l but for only galaxies of high surface brightness with fj, Si g < 26.5 in the models as 
well as the data. 



(2000) , in which they improved the reds hift estimation over 
Fernandez-Soto. Lanzetta & Yahil ( 1999). The difference be- 
tween the results of CSF and DSFO is very small for bright 
galaxies with I&m < 26 (panels a and b). While the CSF model 
reproduces the observed redshift distribution over an en- 
tire range of apparent 7si4-magnitudes considered, the DSFO 
model predicts too many high-redshift galaxies for 4 14 > 26 
(panel c), especially at z > 3. This is because the shorter 
SF timescale in the DSFO model induces earlier formation 
of galaxies compared with CSF. These results are consistent 
with our previous results. We show the same model as CSF 
but for (78 = 1 ■ Since the formation epoch of dark halos shifts 
to higher redshift, a little more galaxies are formed at higher 
redshifts. We also find that the DSF2 model reproduces the 
observed redshift distribution with negligible difference from 
the CSF model with a% = 1 . 

5.3. Lsophotal Area-Magnitude Relation 

Figure [22] shows the isophotal area of galaxies plotted 
against their ^'-magnitude, for which the same observational 



condition employed in the SDF survey is used to calculate 
the isophoto in the SF models. The solid line indicates the 
the mean relation with errorbars of la scatter, predicted by 
the CSF model. Results of other SF models are almost the 
same, and are not shown. The data plotted by the crosses are 
those for the SDF galaxies that are detected in both the K 1 - 
and /-bands. As stressed in our previous papers, the selection 
effects from the cosmological dimming of surface brightness 
of galaxies cannot be ignored in the SAM analysis of galaxy 
number counts. This indicates that the size of high-redshift 
galaxies must be modeled properly. 

With the selection effects correctly taken into account, the 
predicted size should converge towards the limiting magni- 
tude, because faint galaxies with larger area have surface 
brightnesses below the detection threshold and then remain 
undetected. We find from this figure that our SAM galax- 
ies well reproduce the observed area-magnitude relation, and 
are consistent with the SDF galaxies, only when the selection 
bias against faint galaxies with high redshift and/or low sur- 
face brightness is taken into account in the analysis. 
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ure^] but for the SF models without the effects of dynamical response to 
starburst-induced gas removal. 
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ies. Same as Figure fT5l but for the SF models without effects of dynamical 
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FIG. 15. — Distribution of B-band mass-to-light ratio for elliptical galaxies 
for the models of (a) CSF, (b) DSF2, (c) DSF1, and (d) DSFO, where the 
effects of dynamical response to starburst-induced gas removal are ta ken into 
account. Contours and symbols have the same meanings as in Figure [HI but 
for the outermost contour level added showing 10~ 3 times the largest number 
of galaxies in grids. 



5.4. Tully-Fisher Relation (TFR) 

Figure [23] shows the /-band TFR. The thick solid line in- 
dicates the predicted TFR with errorbars of la scatter for 
gas-rich spirals in the CSF model, having more than 10% 
mass fraction of cold gas in the total galactic mass, that is, 
M cn ]d I ( M» + M Kn irt) > . 1 . This criterion is the same as that 
used by ICole et alJ d2000t) . The thick dashed line is for all 
spirals. The tendency that gas-rich spirals for a given mag- 
nitude are always brighte r than gas-poor spirals has already 
been shown in lCole et alJ ( I2000I) . Results of other SF models 
are not shown because of negligible difference. 

The thin dashed, dot-dashed and dotted lines 
are the best-fit results to t he o bserv ed TFRs by 
Math ewson. Ford & Buchhornl (fT99l . iPierce & T ullv 



( 119921) . and lGiovanelli et aT] ( 119971) . respectively. We assume 
that the line-width W is simply twice the disk rotation 
velocity V 10t as usual. For mo del g alaxies, we set Vd to be 
equal to V m t, as mentioned in § 32.51 

The predicted TFR slope for galaxies with small rotation 
velocity V rot < 80km is steeper than that observed, while the 
predicted slope and magnitude for galaxies with larger veloc- 
ity agree well with the observations. The TFR slope is deter- 
mined by ahot, because this parameter relates the mass frac- 
tion of stars in each galaxy, which gives luminosity, to the 
rotation velocity. We use = 4 throughout as a standard, 
but smaller ahot gives shallower slope. For the purpose of 
comparison, we also plot the CSF model with a^ot = 2 by the 
dot-dashed line. In t his case the slope becomes closer to that 
observed. Although Col e et all (|2£)00) adopted ahot = 2 and 
claimed agreement of their model with the observed TFRs, we 
found that this model predicts too many high-redshift galaxies 
to be reconciled with the observed number counts. Moreover, 
mean metallicity of stars in dwarf spheroidals is too high to 
b e con sistent with their observed metallicities, as shown in 
§ 34.41 In order to reproduce both observations, it might be 
worth relaxing our assumptions in estimating the disk rotation 
velocity and/or other physical processes such as SN feedback. 
For example, the SN feedback parameter (3, which is assumed 
to be constant, might depend on whether star formation is ei- 
ther continuous or burst-like. Furthermore, it might evolve 
with redshift. Dynamical response to gas removal, of course, 
provides a promising effect on disk rotation velocity. These 
possibilities should therefore be investigated in more detail. 

5.5. Color-Magnitude Relation ( CMR) 

Figurel24l shows the V — K color versus magnitude relation 
of cluster elliptical galaxies embedded in extended dark halo 
with Vdrc = 10 3 km s . The thick solid, dashed, dot-dashed, 
and dotted lines represent the predicted CMRs in the SF mod- 
els of CSF, DSF2, DSF1, and DSFO, respectively. The CMR 
for each model is obtained by averaging 50 realizations. Er- 
rorbars to the CMR denote the la uncertainties in the DSF2 
model, which are comparable to CSF and DSF1 but is slightly 
larger than DSFO. The thin dashed line is the observed CMR 
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FIG . 17. — Distribution of B-band mass-to-light ratio for elliptical galaxies. Same as Figure [T5l but for only galaxies of high surface brightness with fi e g < 26.5 
in the models as well as the data. 



for ga laxies in the Coma cluster by iBower. Lucev & Ellisl 
(1992) and the thin solid li ne is the same but for the aperture- 
corrected CMR bv lKodama et alJ d!998h . 

All the SF models are not reconciled with the observed 
CMR, although the colors at My - 5 log h ~ -20 are con- 
sistent with the observation. This is a consequence of 
adopting the SN feedback with (ahot,Vta) = (4, 180km s -1 ), 
tog ether with the satellit e -satell ite m erger. As discussed 
bv iKauffmann & Chariot! ( 119981) and lNagashi m7& Goudal 
(2001), a combination of smaller ahot(— 2) and larger Vh ot is 
known to give a better fit to the observed CMR. 

Like the TFR, the faint-end slope of the CMR is deter- 
mined b y cthot, because the CMR is primarily a metallicity se- 
quence ( Kodama & Arimoto and chemical enrichment 
of dwarf galaxies is determined by ahot- The bright-end slope 
of the CMR is nearly fiat, because the SN feedback becomes 
negligible. If we adopt a large value of Vh ot ~ 280 km s" 1 , this 
flat region moves brightwards. Thereby, an expected slope 
from ahot is realized over an almost entire range of magni- 



tudes. Additionally, the satellite-satellite merger makes dwarf 
spheroidal galaxies brighter, while keeping their color un- 
changed. This shifts the CMR bluewards when seen at faint 
magnitudes. 

With these considerations, we can improve the fit to the ob- 
served CMR by adjusting the SN feedback-related parameters 
(othau Vhot), and the merger-related parameter f^. However, 
appropriate choice of their values that could explain the ob- 
served CMR seems to invalidate other successes of our SAM. 
This difficulty is clearly the area of future investigation. 

Furthermore, we find that a bright portion of the CMR is 
also affected by the cooling cutoff V cut . Figure [25] shows 
the V cut -dependence of CMR. The solid and dashed lines rep- 
resent the CSF model with V cut = 250 and 300 km s" 1 , re- 
spectively. The dot-dashed and dotted lines are the same but 
for the DSF0 model. If we adopt larger V cut , the metallicity 
becomes larger because chemical enrichment continues until 
lower redshift. Note that the difference between the models 
of CSF and DSF0 is caused mainly by the difference of mean 
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FIG. 18. — Distribution of mean stellar metallicity for elliptical galaxies in the models of (a) CSF, (b) DSF2, (c) DSF1, and (d) DSFO, where the effects of 
dynamical response to starburst-induced gas removal are taken into account. Contours and symbols have the same meanings as in Figure|9| Only galaxies of high 
surface brightness with /i c s < 26.5 a re shown, be cause the theoretical results in this figure are not sensitive to this threshold . The logarithmic iron abundance 
[Fe/H] is used for the data compiled by Mateo 1 1998) and those translated from Mg2 index by Bender, Burstein & Faber 1 1993), whereas the L^-weighted average 
of logarithmic metal abundances of stars (log(Z„ /Zq))l b for model galaxies. 



stellar age, that is, stars in the DSFO model were born earlier 
and thus older than CSF. Since the CMR is sensitive to Vc Ut , 
more knowledge of the gas cooling is needed to fix the cutoff 
on cluster scales. 

This situation might be improved by taking the following 
procedures. First is to use a lower value of erg below unity 
as rec ent observations suggest ers — 0.8 (e.g., Soer gel et alJ 
2003). This value provides a statistically late epoch for den- 
sity fluctuations to collapse. Therefore, a large value of V cut 
does not result in the formation of monster galaxies. Second is 
the aperture effects. Kawata & Gibson (2003) showed the im- 
portance of aperture effects by using their chemo-dynamical 
simulations. The predicted V—K colors of their giant elliptical 
galaxies are on the average 3.2 at My -5 log h ~ -22 and agree 
well with the observed colors within 5 kpc aperture. How- 
ever, the observed V—K color within 99kpc aperture is nearly 
equal to 3.0, which is even bluer than our predicted color at 



the same m agnitude. Third is the effects of UV background. 
Nagashima & Gouda (2001) showed that the photoionization 
by the UV background has a similar effect to the SN feedback. 
Therefore, introduction of the photoionization is equivalent to 
adopting larger Vh ot . Further investigation on the CMR as well 
as the TFR are needed alo ng these lines. 

It should be noted that Okamoto & Nagashima ( 2003) have 
successfully reproduced the CMR slope by using a SAM com- 
bined with an A^-body simulation. Their SF model corre- 
sponds to CSF with ahot = 2 and Vh ot = 200 km s . As ex- 
plained above, however, adopting this small value of «hot fails 
to explain other observations unless the merger strength is ad- 
justed. In this sense, high resolution Af-body simulations is 
highly awaited to follow a full trace of merging histories of 
progenitor halos of clusters. 

6. COSMIC STAR FORMATION HISTORY 
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FIG. 19. — Galaxy number counts in the (a) U300, (b) B450, (c) Vgo6 and (d) /814 bands, respectively. The thick and thin lines represent the SF models with and 
without the selection effects, respectively, arising from the same detection threshold of surface brightness of galaxies as employed in the HDF survey. For each 
of the cases, the solid and dashed lines represent the models of CSF and DSF2, respectively, where the effects of dynamical response to starburst-induced gas 
removal are taken into account. Open circles with errorbars indicate the observed HDF counts. Other symbols are ground-based observed counts. 



Recently the cosmic SF history, which is a plot of SFR in a 
comoving volume agai nst redshift, is wide ly used to examine 
the global SF history (jMadau et alJU996l) and also in SAM 
analy ses dBaugh et alJ 119981 iSomerville. Primack & FabeJ 
12001 . Since the SFR is very sensitive to the SF timescale, 
in Figure[26]we plot the redshift evolution of cosmic SF rate 
for the four models with different SF timescales. Thick lines 
denote the total SFR and thin lines the SFR only for star- 
burst. Symb ols with errorbars indi cate observational SFRs 
compiled by Ascasibar et al. (2002). While there are large 
scatters between individual data points, the CSF and DSF2 
models broadly agree with the observations. 

As expected from the redshift dependence of SF timescale, 
the CSF model predicts a maximum SFR at relatively low red- 
shift, z < 2, and the redshift at which this maximum occurs 
becomes larger, in order from DSF2 (dashed line) via DSF1 
(dot-dashed line) to DSF0 (dotted line) as the SF timescale 
becomes shorter at high redshift (see Figure|2j. 



The fraction of cold gas available for starburst becomes 
larger for longer SF timescale at high redshift. Because of 
the same reason as above, the CSF model by the thin solid 
line gives the largest fraction of cold gas at major merger and 
therefore the highest SFR for starburst among the four models 
under consideration. Accordingly, the dynamical response is 
the most effective in the CSF model. 

7. COOLING DIAGRAM REVISITED 

ISilkl ( 119771) and iRees & OstrikeJ (119771) proposed the so- 
called cooling diagram, that is, the distribution of galaxies on 
the density versus temperature diagram. Since then, a num- 
ber of authors used this diagram as crucial constraints on the 
formation of galaxies i n the framework of monolithic clou d 
collapse scenario (e.g., Faber 1982; Blumenthal et al. 1984). 
For example, characteristic mass of galaxies can be evaluated, 
because the evolutionary path of gas clouds in this diagram 
gives an initial condition of density and temperature expected 
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Fl G. 20. — Galaxy number counts in the K' (2. 13fj,m) band. Same as Figure 
fT9l but for the near-infrared counts in the SDF. In addition to the models of 
CSF (solid lines) and DSF2 (dashed lines), the results of the DSF0 model 
(dot-dashed lines) are shown for comparison. 
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FIG. 21.— Redshift distribution of the HDF galaxies for (a) 22 </ 814 < 24, 
(b) 24 < 7 8I4 < 26, and (c) 26 < 7 8 i 4 < 28. The thick solid, dashed and 
dot-dashed lines represent the models of CSF, CSF with a% — 1 and DSF0, 
respectively, where the effects of dynamical response to starburst-induced gas 
removal are taken into account. The histogr am in each pa nel is the observed 
photometric redshift distribution by Furusawa et al. 1 2000). 
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FIG. 22. — Isophotal area of the SDF galaxies against ^'-magnitude. 
The solid line indicates the the mean relation with errorbars of lcr scatter, 
predicted by the CSF model, where the effects of dynamical response to 
starburst-induced gas removal are taken into account. The isophoto in this 
model is calculated from the observational condition employed in the SDF 
survey. The crosses indicate the observational data from the SDF survey. 
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FIG. 23. — /-band Tully-Fisher relation of spiral galaxies. Variants of 
the CSF model with the effects of dynamical response are shown for the 
cases of only gas rich spirals having the cold gas more than 10% in mass 
(thick solid line), all spirals with our standard set of parameters (thick dashed 
line), and gas rich spirals with a smaller value of a^ot = 2 (thick dot-dashed 
line). Errorbars indicate lcr scatter around the predicted mean. The observed 
TFRs are shown as a mean relation given by Mathewson. Ford & Buchhorn 
(1992T) fthin dashed lin e), IPierce & Tufhl 119921) (thin dot-dashed line), and 
IGiovanelli et alj 11991) (thin dotted line). 



from density fluctuation spectrum in the early universe. 

Figure^] shows the contours of galaxy distribution in the 
cooling diagram for ellipticals (top panel), lenticulars (middle 
panel), and spirals (bottom panel). The velocity on the hori- 
zontal axis indicates the velocity dispersion of bulge compo- 
nent for elliptical and lenticular galaxies, and the disk rotation 
velocity for spiral galaxies. The baryon density on the vertical 
axis is estimated by dividing the baryonic mass in individual 
galaxies by their volume. We simply assume a sphere of effec- 
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FIG. 24. — V -K color-magnitude relation of elliptical galaxies in clusters. 
Such galaxies are sampled from the clusters having V c = 10 3 km s , and 
the theoretical CMRs are shown for the models of CSF (thick solid line), 
DSF2 (thick dashed line), DSF1 (thick dot-dashed line), and DSFO (thick 
dotted line), where the effects of dynamical response to starburst-induced gas 
removal are taken into account. The observed CMR for the Coma cluster is 
shown as a mean relation by Bower, Lucey & Ellis 1 1992) (thin dashed line), 
and the aperture-corrected CMR for the same cluster by Kodama et al. 1 1998) 
(thin solid line). 
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FIG. 25. — V — K color-magnitude relation of elliptical galaxies in clusters. 
Same as Figure l24l but for the effect of changing V cu t. The thick solid and 
dashed lines show the theoretical CMRs in the CSF model with V cllt = 250 
(our standard choice) and 300 km s _1 , respectively. The dot-dashed and dot- 
ted lines similarly show these two cases for the DSFO model, respectively. 
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FIG. 26. — Cosmic Star Formation Histories. The solid, dashed, dot-dashed 
and dotted lines denote CSF, DSF2, DSF1 and DSFO, respectively. The thick 
lines indicate the total SFR and the thin lines the SFR only for starburst. 
Symbols with errorbars are observational data compiled by Ascas lbar et alj 
(2002). 



solar compositions, respectively. These timescales are 

kT 



<r - 3 Pg 
'cool — 



2//Wp7i|A(7 , ,Z gas )' 
= (24^G Aot r I/2 , 



(21) 
(22) 



where n e is the electron number density, p gas is the gas den- 
sity, p to t is the total mass density including dark matter, and 
A(r,Z gas ) is the cooling function depending on both the gas 
temperature T and metallicity Z gas . Wherever necessary be- 
low, we transform the gas temperature into velocity dispersion 
based on the virial theorem, and identify the gas density with 
baryon density. The thin solid curves represent the sequences 
of galaxies originated from the CDM density fluctuations with 
density contrast of la and 3a, respectively. It is evident from 
this figure that our SAM galaxies are distributed within a re- 
gion of r coo i < T grav and their morphologies are distinctly seg- 
regated from each other, as observed in this diagram. 

Figure shows the distribution for all galaxies without 
imposing any selection bias against low surface brightness. 
Comparison of Figures l27l and 1271 indicates that many dwarf 
galaxies have low surface brightness and are distributed to- 
wards a region of r coo i > T glav characterized by low baryon 
density and low circular velocity. This extended distribution is 
prominent only for elliptical galaxies for which the starburst- 
induced gas removal followed by dynamical response has 
the maximal effect. Such galaxies of low surface brightness 
would have been detected in recent ultradeep surveys where 
the detection threshold is set below = 26.5. 



tive radius for elliptical and lenticular galaxies, and a cylinder 
of scale height being one tenth of the effective radius for spi- 
ral galaxies. The contours are plotted excluding low surface 
brightness galaxies with fi Ci B > 26.5. 

The dashed and dot-dashed lines represent the cooling 
curves that the cooling timescale is equal to the gravitational 
free-fall timescale r coo i = T grav for the cases of primordial and 



8. SUMMARY AND CONCLUSION 

We have investigated the formation and evolution of galax- 
ies in the context of the hierarchical clustering scenario by us- 
ing the Mitaka model, or our SAM in which the effects of dy- 
namical response on size and velocity dispersion of galaxies 
are explicitly taken into ac count, according to the formula by 
Naaashima & Yoshii J2003I) for galaxies consisting of baryon 
and dark matter. This paper is therefore an extension of previ- 
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FIG. 27. — Cooling diagram for (a) elliptical galaxies, (b) lenticular galax- 
ies, and (c) spiral galaxies. Contours indicate the theoretical distribution pre- 
dicted by the CSF model with surface brightness cutoff at fi e g = 26.5 and the 
effects of dynamical response taken into account in the analysis. The levels of 
contours are the same as in Figure^] Top panel shows a broad distribution of 
elliptical galaxies consisting of dwarf ellipticals on the left and normal/giant 
ellipticals on the right, which are bifurcated reflecting a variety of star forma- 
tion and merger histories. The dashed line represents r coo i = T grav for the gas 
of primordial chemical composition, and the dot-dashed line for the gas of so- 
lar chemical composition. The thin lines represent two sequences of collapse 
of overdense regions with la and 3a fluctuations against CDM spectrum. 



ous a nalyses in the context of monolithi c cloud co llapse sce- 
nario (iDekel & Silkll986tlYosh"ii & Arimotoli987i) . 

A A-dominated flat universe, which is recently recognized 
as a standard, is exclusively used here. The investigation 
mainly focuses on elliptical galaxies which are assumed to 
be formed by major merger and starburst-induced gas removal 
followed by dynamical response of the systems. While a mass 
fraction of removed cold gas, after heated, is determined by 
the circular velocity of galaxies similar to the traditional col- 
lapse models, the total amount of cold gas depends on the 
formation history of galaxies which is realized by the Monte 
Carlo method based on the power spectrum of density fluctu- 
ation predicted by the CDM model. 

The Mitaka model, which we have constructed in this paper, 
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FIG. 28. — Cooling diagram. Same as Figure l27l but for all galaxies without 
imposing any selection bias against low surface brightness. 



is found to reproduce a wide variety of observed characteris- 
tics of galaxies, particularly their scaling relations among var- 
ious observables such as magnitude, surface brightness, size, 
velocity dispersion, mass-to-light ratio and metallicity. This 
strongly supports the CDM cosmology and merger hypothesis 
of elliptical galaxy formation even on scales of dwarf galax- 
ies. 

Most of model parameters related to star formation, SN 
feedback and galaxy merger are constrained by the local lumi- 
nosity function and mass fraction of cold gas in spiral galax- 
ies. As an extra parameter to be furthermore constrained, we 
have examined redshift dependence of SF timescale through 
comparison among the four SF models denoted by CSF, 
DSFO, DSF1, and DSF2. The CSF refers to a constant 
SF timescale against redshift and has been suggested to 
be consistent wit h observations of galaxy number counts 
(Nagashimaetal. 2001, 2002), quasar luminosity function 
( Kauffmann & Haehnelt 2000| lEnoki. N agashi ma & Goudal 
2003), and number and metallicity evolution of damped Ly- 
a sys tems (ISomerville. Primack & Faberll2001t lOkoshi et alJ 
2002). The DSFO refers to a SF timescale simply proportional 



22 



NAGASHIMA & YOSHII 



to dynamical timescale. The DSF1 and DSF2 are intermedi- 
ate between the CSF and DSFO (see Figure |2j. It is found 
that the DSFO model fails to reproduce observed properties of 
local dwarf spheroidals and galaxy counts. Among the above 
four SF models, the CSF model has the longest SF timescale 
at high redshift and therefore the largest amount of cold gas at 
that epoch. Accordingly, the fraction of removed gas during 
starburst is the largest in the CSF model, giving many ex- 
tended galaxies in agreement with observation. In this sense, 
the model of constant star formation (CFS) or at most mildly 
evolving star formation (DSF2) is favorable. Note that dom- 
inant dark halo, in which size and velocity dis persion do not 
change during gas removal as considered by iDekel & Silk 
( 1986), apparently explains the observed scatter of size. How- 
ever, arguments from physical ground indicate that the dy- 
namical response should play a significant role in forming 
galaxies with low velocity dispersion of o-q(ID) < 10km s -1 , 
as actually observed in the Local Group. 

There are some areas of further improvements in the current 
Mitaka model. Because of strong dependence of SN feed- 
back on circular velocity assuming ahot = 4, dwarf galaxies in 
our SF models are too faint to agree with expected magnitude 
from observed TFRs. Furthermore, because of the same rea- 
son, dwarf galaxies are too blue to agree with expected colors 
from observed CMRs for cluster elliptical galaxies. On the 
other hand, such large value of «hot stops chemical and pho- 
tometric evolution of elliptical galaxies at early epochs and 
is required to explain their observed low stellar metallicity as 
well as the galaxy number counts in the HDF and SDF In or- 
der to settle these contradictions, some new ingredients need 
to be introduced in the SAM analysis. At least, for example, 



we must know how efficient the gas cooling is in massive dark 
halos and whether SN feedback at starburst works in the same 
way as that in disks. 

We adopted a fitting mass fun ction of dark halos by 
Yahagi. Nagash ima & Ybshiil J2003h for given redshift instead 
of the often used PS mass function. This YNY mass func tion 
is slightly different from those by Sheth & Tormenl ( 119991) and 
iJenkins et alJ d2001l> . but is conf irmed to provide a better fit to 
recent Af-body results given by Yahagi & Yoshii ( 200l|) and 
Yahagi (2002), which is discussed in more det ail in a sep- 
arate paper ( Yahagi, Nagashima & Yoshii 2003). Since the 
number density of dark halos affects SN feedback-related pa- 
rameters to be chosen, it is very important to establish the 
mass function of dark halos, particularly on small mass scales, 
by high resolution iV-body simulation. Although such SAM 
analyses have just begun (e.g., Ben son et alJ200ltlHellv et al.l 
2003a b), further investigation will obviously be required. We 
are upgrading the Mitaka model accommodated with full high 
resolution N-body simulations, and the results will be given 
elsewhere in the near future. 
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APPENDIX 

DYNAMICAL RESPONSE TO GAS REMOVAL 

We show the dynamical response of size and velocity dispersion to supernova-induced gas re moval in the two-component 
galaxies consisting of baryon and dark matter. The details are given in Nagashima & Yoshii (2003). 
All elliptical galaxies in this paper have a density profile similar to the Jaffe model ( Jaffe 1983), 



p(r)-. 



r 2 (r+rt,) 2 ' 



(Al) 



where pi, and ri, are the characteristic density and radius, respectively. This profile well approximate the de Vaucouleurs r 1 / 4 
profile for stars. The effective (half-lig ht) radius r e defined in the projected surface is related to the half-mass radius rt in three 
dimensional space as r e = O.lAArb (Nagashima & Yoshii 2003). For the distributi on of dark matter, the singular i sothermal distri- 
bution (pd oc r" 2 ) is considered. We found that the Navarro-Frenk- White profile (iNavarro. Frenk & White! 19 97) provides almost 
the same results as those given by singular isothermal sphere. Note that according to recent high resolution iV-body simula- 
tions, which have reveals that many substructures survive even in v irialized halos, we assume that subhalos exist as underlying 
gravitational potential for satellite galaxies. As mentioned in § 32.21 these subhalos are tidally truncated. 

Because the derivation is long and complicated, we describe it only briefly. Just after the merger, the system is virialized 
immediately. Then the starburst occurs and a part of gas is gradually removed. The size and velocity dispersion change following 
the gas removal adiabatically. During the gas removal, the dark matter distribution is assumed to be not affected. 

In the fa llowings we d erive useful simplified formulae for the dynamical response. As shown in Appendix B of 
Nagashima & Yoshii] d2003h . the relationship between density and size for the adiabatic gas removal is given by 



Y = 



h(yi,Zi)p(Zf) + q(Zf) 



(A2) 



where Y = y/fyi, = Pi/pj, Z;/ = r^f/r^, and the function h,p and q are 



h (y, z) = yz 4 + - [-z + ( 1 + z) ln( 1 + z) ] , 



P(z) = 



z 4 ' 



q(z) = -^[-Z + (l+z)m+z)]. 



(A3) 
(A4) 
(A5) 
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FIG. 29. — Dynamical response of galaxy size to removed mass from the 
system. The approximate results given by equation (A9) are shown for D = 
0.05 and 0.2 ( solid lines). For each value of D, the exact results given by 
equation IA6I are shown for zi = 0.5 (dashed lines) and 0.05 (dot-dashed 
lines). 



In order to obtain the size after gas removal, we need to know the inverted function of the above and it will be very complicated. 
Therefore, we approximate the above function. The equation can be reduced to the following expression, 

r i (A6) 

R 4 2y iZ 4 R 4 

where R = Zf/zi- Expanding t/; so (/?,z,) around R = 1 and z, = and picking out the lowest order term, we obtain 

Y=\~- R -4, (A7) 

R 4 2 R 4 ' ' 

where D = 1 /yizf- Defining the ratio of final to initial masses, A4 = Mj/Mj = YR 3 , the above equation is transformed to 

A4 = ±-? — ■ (A8) 



Then, inverting this, we obtain 



R 2 R 
l+D/2 

R -M^f2- (A9) 



By using the virial theorem, the velocity dispersion is 

o_ L= l YR 2 + Df(z f )/2 
at V l+0/fe)/2 ' 

where 



(A10) 



/(z)= l"£±£) +ln f 1 + I, (All) 



z \ z, 

Figure [29] shows the relation between M. and R based on this approximate formula with D = 0.05 and 2 (solid lines). For 
each value of D, the exact solution is also shown for two cases of z, = 0.05 (dot-dashed line) and 0.5 (dashed line). Clearly our 
approximate relation follows the exact solution quite well. 
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